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WAVE  RADIATION  FROM  SLENDER  BODIES 


J.  N.  Newman 

DEPARTMENT  OF  OCEAN  ENGINEFRING 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
CAMBRIDGE,  MASS. 


ABSTRACT 

Wave  radiation  by  a slender  body  is-  discussed,  with 
emphasis  on  the  scale  of  the  characteristic  wavelenath 
in  relation  to  the  disparate  width  and  lenath  of  the 
body.  The  discussion  is  illustrated  by  reviewina  the 
analysis  of  oscillatory  forced  motions  of  a slender 
floating  body  in  otherwise  calm  water.  Two  ccemp- 
lementary  regimes  are  considered,  where  the  wavelength 
is  long  (comparable  to  the  body  length)  or  short 
(comparable  to  the  beam)  . The  limitations  o.-f  these  two 
theories  suggest  the  need  for  a unified  approach,  where 
geometric  slenderness  is  assumed  without  restricting 
the  wavelength.  A possible  aonroach  to  the  uivified 
theory  is  outlined  in  general  terms,  and  illustrated 
in  detail  by  considering  the  simpler  ohysical  oroblem 
of  axisymmetric  acoustic  radiation  bv  a slender  body 
in  an  unbounded  medium. 
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Feinier  Tirntran  influenced  zka  field  of 
applied  r<atke~.atias  by  kis  ovn  vovk,  that 
of  kis  students  at  Felft,  and  tkrcuok  t'ne 
•jork  of  rally  others  who  uere  kis  students 
in  a generalized  sense.  The  osoasional 
opportunities  to  discuss  -w  o'on  •ocr<  ■citk 
Prefessor  Ti-ran  were  invaluable.  Traces 
cf  kis  s t i~:< a t i ny  i'ls'ykr  arc  rrrs-nt  in 

Py-:  jsi:,~,  wk'ok  serves  bc'.k  ' : krn.'ur 
and  to  yerre  ; :,a  te  kis  u'rltue  i flu  . 


1 . INTRODfCmiO:: 

The  field  o^  slencer-shio  thcorv  is  a 
loaical  develocrent  ‘'ron  the  earlier 
analysis  ot  slender  bodies  in  aerodynamics, 
ftov.-eyer  a sioni-^icant  comolication  results 
for  shins  and  ether  'loatino  slender 
bodies,  due  to  the  extra  lenoth  scale 
assoriat-oJ  v,-ith  the  characteristic  wave- 
lonnth  (■•  !,  ~h;' •?!■,  fast  be  related  to  the 
disrarate  Icnmtl,  (L)  and  beam  (B)  o*  the 
vosse)  . ^'n  analc'.-ous  situation  exists  <'or 
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Acoustic  radiation  in  an  unbounded  r.iodium. 
With  the  fundamental  geometric  slender- 
ness parameter  B/L=e<<l,  we  seek  an 
asymptotic  solution  valid  to  leading 
order  in  e . For  this  purpose  it 
generally  is  necessary  to  restrict  the 
order  of  magnitude  of  the  wavelenqth  to 
one  of  two  complementary  regimes.  For 
long  waves,  comparable  to  the  body 
length,  X/L=0(1),  whereas  in  the  short- 
wavelength  case  X/B=0(1).  (Some  workers 
reserve  the  term  "slender-body  theory"  or 
"slender-ship  theory"  for  the  long-wave- 
length regime.  A more  liberal  definition 
is  adopted  here,  the  word  "slender"  being 
applied  only  to  the  geometry  of  the  body.) 

In  the  long-wavelength  regime,  inter- 
actions are  significant  between  adjacent 
sections  of  the  body,  but  wave  effects 
are  absent  from  the  near  field  close  to 
the  body  surface.  For  the  ship-motion 
problem,  the  importance  of  transverse 
gradients  in  the  inner  region  is  such 
that  the  linear  f ree-surface  boundary 
condition  degenerates  to  the  "rigid-wall" 
condition  of  zero  vertical  velocity. 

This  greatly  simplifies  the  theory,  but 
also  restricts  its  domain  of  applicabil- 
ity to  sub-resonant  frequencies. 

In  the  short-wavelength  regime  wave 
effects  are  present  in  the  near  field, 
but  one  recovers  a simple  strip  theory 
without  interactions  between  sections. 

In  one  sense  this  result  is  very  satisfy- 
ing, as  it  provides  some  rigor  to  the 
otherwise  empirical  strip  theory  of  ship 
motions.  However,  the  consistent  leading- 
order  theory  is  again  oversimplified  for 
practical  purposes.  In  particular,  the 
hydrodynamic  effects  of  the  ship’s  for- 
ward speed  are  absent  completely  from  the 
leading-order  theory,  as  suqgested  by 
Vossers  (1962)  and  confirmed  by  Joosen 
[1964]  and  Cgilvie  and  Tuck  (1969). 

Faced  with  the  choice  between  these  two 
complementary  theories,  and  the  obser- 
vation that  practical  ship  motions  occur 
in  wavelengths  which  generally  are  inter- 
mediate between  the  length  and  beam. 


it  is  natural  to  seek  a "unified"  slender- 
body  theory  where  the  wavelenoth  is  un- 
specified in  advance.  A suggested  approach 
to  that  objective  is  tne  principal  con- 
tribution of  this  lecture. 

First,  I shall  outline  briefly  the 
separate  treatments  of  lone-  and  short- 
wavelenoths  for  the  shio-motion  problem. 

A possible  approach  is  then  developed  for 
the  more  qeneral  case,  and  illustrated  bv 
treatinq  the  analoaotis  problem  of  acoustic 
radiation  from  a slender  bodv  in  an  un- 
bounded medium. 

For  pedagogic  reasons  the  discussion  will 
be  simplified  as  much  as  possible.  A 
significant  restriction  is  to  consider 
only  the  radiation  problem,  with  a slowly- 
varying  distribution  of  the  prescribed 
normal  velocity  on  the  body  surface.  The 
diffraction  problem  for  scattering  of 
incident  waves  can  be  solved  in  a similar 
manner  if  the  wavelencth  is  larce,  but 
for  short  wavelenqths  the  diffraction 
problem  is  more  difficult. 

Another  restriction  is  to  neglect  forward 
speed  in  the  shio-motion  problem.  This 
simplifies  the  boundary  conditions  of  the 
problem,  and  the  resultino  solution.  The 
effects  of  forward  speed  are  included  in 
the  long-  and  short-wavelenath  problems, 
respectively,  by  Newman  and  Tuck  (1964) 
and  Cqilvie  and  Tuck  (1969).  No  fundamen- 
tal difficulties  are  anticipated  in  adcot- 
inq  t.he  suggested  unified  approach  to 
include  forward-speed  effects. 

2.  THF  SHIP-MOTION  PROBLEM 
The  problem  of  interest  here  is  to  analyse 
the  motions  of  a ship,  or  more  oenerallv 
of  any  elonqated  vessel  floating  on  the 
free  surface.  Cartesian  coordinates 
(x,y,z)  are  defined  with  z=0  the  plane  of 
the  undisturbed  free  surface,  and  the 
direction  of  the  positive  z-axis  upward. 

It  is  convenient  to  nondinensionalizs  and 
orient  these  coordinates  such  that  the 
longitudinal  bodv  axis  occupies  the  seg- 
ment (0,1)  of  the  x-axis. 

The  notions  of  the  body  and  surrounding 
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fluid  are  assumed  to  be  oscillatory  in 
time,  with  the  complex  factor 
suppressed.  The  motions  are  assumed  also 
to  be  of  small  amplitude  by  comparison  to 
the  wavelength  and  body  dimensions.  Thus 
a linearized  boundary-value  problem  can  be 
justified,  after  makina  the  usual  assump- 
tion that  the  fluid  is  ideal  and  incompres- 
sible. 

The  fluid  velocity  vector  is  expressed  as 
the  positive  gradient  of  a potential  i (x) , 
which  is  governed  throughout  the  fluid 
domain  by  the  Laplace  equation 


^ as  (J)  + 

^ ^XX 


The  two  principal  boundary  conditions  are 
on  the  wetted  surface  of  the  body  and  on 
the  free  surface.  The  appropriate 
boundary  condition  on  the  body  surface  S 
is 

34/3n  = V , 


where  is  the  specified  normal 

velocity.  In  the  linearized  theory,  (2) 
is  imposed  on  the  mean  position  of  S , 
and  the  free-surface  boundary  condition 
takes  the  form 

KCi-i^  = 0,  on  2 = 0 , 

(3) 

with  K=aj*L/g  . For  deep  water  the  fluid 
motion  should  vanish  as  z , and 

K=2^L/).  is  the  (nondimensional)  wavenumber 
of  the  radiated  waves. 

This  boundary-value  problem  is  comoleted 
by  imposing  the  radiation  condition 

e , as  R=  (x' +y ' ) Vi  oo 


by  a suitable  initial-value  oroblem,  with 
a state  of  rest  at  t=  . A convenient 
aooroach  is  to  retain  the  time-deoendence 
e^^^  , but  with  a comolex  freouencv  u 
having  a vanishingly  small  negative  cart, 

0)  -*■  CD  -iO  . (5) 

From  the  definition  of  the  wavenumber  it 
follows  that 

K ■*  K -iO  . (6) 

We  shall  use  (6)  in  place  of  (4) , since 
this  avoids  the  asymptotic  analysis  to 
verify  the  far-field  wave  form.  Further- 
more, (5)  and  (6)  can  be  utilized  in  the 
acoustic  problem  of  Section  3. 

The  body  geometry  is  characterised  by  a 
length  LHl,  beam  (B)  and  draft  (T) , such 
that  {B,T)  = 0(e)  . Our  aporoach  is  based 
on  the  method  of  matched  asvmpto.tic  ex- 
pensions.  Thus,  we  define  an  inner  region 
adjacent  to  the  body  surface,  where  • 

(y,2)  =»  0(e)  , and  an  outer  region  where 
(y,z)  =0(1)  . By  assumption  there  exists 
an  overlap  region  E<<(y,z)<<l  , where  the 
solutions  obtained  separately  in  the  inner 
and  outer  regions  are  both  valid  and  can 
be  matched.  This  technique  is  described 
by  Van  Dyke  (1975)  , and  was  first  applied 
to  slender  ships  by  Tuck  (1963) . 

In  the  inner  region,  a coordinate-stretch- 
ing argument  can  be  used  to  show  that 
gradients  in  the  transverse  plane 
(3/3y,3/3z)  are  0(e~’)  , and  dominate 
the  longitudinal  derivative  3/3x=0(l). 
Thus  (1)  reduces  to  the  two-dimensional 
Laplace  equation 


t + V = 0 , 

yy  zz  ' 


with  the  error  a factor  1+0 (e')  . 

The  solution  in  the  inner  region  is  of 
the  form 


Alternatively,  the-  radiation  condition  (4) 
can  be  deleted  if  the  s toad /-state 
oscillatory  t Imo-depcndonce  is  replaced 


t(x,y,3)  rr  : (y,3;x)  + f (X)  . 


(8) 
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Hftre  denotes  the  solution  of  a 

two-dimensional  boundary- value  problem, 
and  f (x)  is  a trivial  solution  of  ( 7) 
which  is  determined  ultimately  by  matching 
to  the  outer  solution. 

In  the  outer  reoion  there  is  no  preferred 
direction,  and  all  three  components  of 
the  Gradient  ooerator  are  of  the  same 
order  of  magnitude.  Thus,  the  three- 
dimensional  Laplace  equation  (1)  holds 
in  the  outer  region,  without  simplifica- 
tion. From  geometric  considerations,  the 
components  of  the  unit  normal  vector  n 
directed  into  the  body  surface  are 


''x  “ ' 


(ny,n^)  =»  0(1) 


(9) 

(10) 


The  estimates  (9-10)  must  be  modified  near 
the  bow  and  stern  unless  these  ends  are 
sharply  pointed.  Hereafter  we  restrict 
the  body  geometry  such  that  (9-10)  ace 
uniformly  valid  along  the  lencth  of  the 
‘ship.  Furthermore,  a longitudinal  source 
distribution  will  be  utilized  subsequently 
to  derive  the  outer  solution,  and  it  will 
be  necessary  to  restrict  the  source 
strength  near  the  ends.  Ke  shall  assume 
that  the  source  strength  is  zero  at  the 
body  ends,  especially  in  connection  with 
various  partial  integrations;  a stronger 
condition  may  be  necessary,  particularly 
if  the  solution  is  to  be  uniformly  valid 
near  the  ends.  This  question  is  dis- 
cussed by  Ursell  (1962) , but  is  not  of 
practical  importance  in  the  leading-order 
analysis  of  vertical  ship  motions. 

It  is  obvious  that  the  body  boundary 
condition  (2)  is  apolied  in  the  inner 
region,  and  the  radiation  condition  (4) 
is  valid  only  in  the  outer  region.  The 
free-surface  condition  (3)  must  be  aoolied 
in  both  the  inner  and  outer  reoions,  and 
the  manner  in  which  this  is  done  daoends 
cn  the  order  of  macnitude  of  the  wave- 
ler.ath  or  wavenumber  K . 


Lono  V>'avG lengths  , K=0(1) 

If  K=0(1),  the  dominance  of  transverse 
gradients  in  the  inner  region  leads  to 
the  "rigid"  free-surface  condition 


on  z=n 


(11) 


in  Place  o^  (3).  Bv  reflection  about  the 

plane  z=0,  t is  the  solution  of  a 

two-dimensional  flow  exterior  to  the 

wetted  contour  C of  the  ship  and  its 

image  C above  z=0  , with  prescribed 

normal  velocity  V on  C and  the  same 
n 

normal  velocity  on  C . 

From  conservation  of  mass  will 

include  a source  strength  proportional 
to  the  net  flux  associated  with  the  normal 
velocity  on  C . Thus  for  heave  and 

pitch 


2i-  (S')  r , 


r=(v*+z*))^2  >>e 


(12) 


The  source  strength  !7(x)  is  the  product 
of  the  local  half-beam  and  vertical  velocity. 

The  remaining  motions  o^  the  ship  will  not 
result  in  a net  source  strength  and 


(2D); 


u (x) 'V  log  r,  r>>c 


(13) 


For  sway  and  yaw  the  dipole  moment  j£(x) 
can  be  related  to  the  added  mass  of  C+C  . 

An  arbitrary  "constant"  f(x)  can  be  added 
to  , and  this  constant  contributes 

to  the  pressure  in  the  inner  reoion  adiacent 
to  the  body  surface.  This  affects  the 
vertical  force  on  the  profile  C,  but  not 
the  horizontal  force.  This  constant  is 
significant  for  pitch  and  heave,  and  must 
be  determined  bv  matching  (12)  with  the 
outer  solution.  By  comparison,  for  swav 
and  vaw,  the  constant  can  be  neglected 
and  the  local  solution  in  the  inner  reoion 
is  given  sinolv  bv  the  two-dimensional 
potential  of  the  form  (13). 
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Hereafter  we  consider  only  pitch  and 
heave.  In  the  outer  region  a distribution 
of  three-dimensional  wave  sources  is 
required,  governed  by  the  complete  free- 
surface  condition  (3),  The  necessary 
source  potential  or  Green  function  is 
well  known,  cf.  Wehausen  and  Laitone  (1960, 
eq.  13.  17).  For  a source  of  complex 
time-dependence  , situated  at  the 

point  (C,0,0)  , the  velocity  potential 
can  be  expressed  in  the  form  G(x-C,y,z), 
where 


tematic  manner  usino  Fourier  transforms, 
followina  Ursell  (1962) . With  an  asterisk 
(*)  used  to  denote  the  Fourier  transform 
with  resoect  to  x,  we  define 

())*(y,z;k)=re^’'^  « (x,Y,z)  dx  , (17) 

/ 

and  similarly  for  other  functions  of  x . 
By  the  convoLation  theorem 

'J'*  = -hr  G*  , (18) 


Glx.y.zl  . 

J^^k' (x^+yM . (14) 

Note  that  (14)  satisfies  the  (three- 
dimensional)  Laplace  equation  (1),  and 
free-surface  condition  (3) . In  view  of 
(6)  the  appropriate  contour  of  inte- 
gration is  above  the  pole  at  k'=K  . 

With  this  choice  of  contour  (14)  can  be 
shown  by  asymptotic  expansion  to  satisfy 
the  radiation  condition  (4).  Near  the 
source  point,'  (14)  is  dominated  by  the 
singularity 

G 5 I”  dk'  ^k'  (x^+y=')’^2je’^'^ 


where  G*  is  the  Fourier  transform  of 
(14)  . 

After  substitution  of  the  intearal  reo- 
resentation  for  , the  Fourier  trans- 

form of  (14)  is  aiven  in  the  form 


G*=i_r 


2T 

0 


gk'z  + ik-(x  COS9+V  Sin6)^g  ^ 


(19) 

The  x-integral  in  (19)  yields  a delta 
function.  After  an  apolication  of  General- 
ized harmonic  analysis  it  follows  that 


. r exp[z(k^+  +i<vl 

^ -”  (k^  + <-  ) Va  -K 


dx 


(20) 


where 


In  the  outer  region  far  from  the  surface 
of  the  body,  the  potential  can  be  ex- 
pressed as  a source  distribution  on  the 
body  axis  with  unknown  strenqth  c(x),  in 
the  form 

i>=  -j:r  a(5)  G(x-c,y,z)  dC  . (16) 

' o 

In  our  subsequent  analysis  it  will  be 
assumed  that  a is  a regular  function, 
which  vanishes  at  the  body  ends. 

The  inner  expansion  cf  (16),  for  small 
(y,n),  is  required  for  matching  with  the 
inner  solution.  This  noarfield  approxi- 
mation of  (16)  can  be  derived  in  a sys- 


< = k'  sin  S 


An  expansion  of  (20)  in  modified  Bessel 
functions  o^  aroument  (Ik'r)  is  derived 
bv  Ursell  (1962)  . ’^or  small  values  oF 
this  arnument  the  leadino-order  terns  are 


G*  (r  sin  e , 


-2y 


+ ^ 

( V-K- 


-r  cos  ®;k)=-2  loa^i!kIr| 

/o-cos”’ (K/' k! ) 1 

^ y.  I -iT-co3h~'  (K/' k I ) / 


for 


(21) 
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(Ursell's  eypression  is  conjugate  to  (21), 
due  to  the  ooposite  convention  for  the 
con'plex  time  dependence,  and  includes 
hiqher-order  terms  proportional  to  Kr.) 

The  inverse  transform  of  (21)  can  be 
evaluated,  following  Ursell  (1962).  After 
a second  application  of  the  convolution 
theorem,  the  inner  expansion  of  the  outer 
solution  can  be  expressed  in  the  alterna- 
tive forms 

^ = a (x)  log  r 

+ T7  [’  a’ (5)L  (x-C)  dC 

“ ' d 


of  Newman  (1964)  . The  equivalence  of  (22a) 
and  (22b)  follows  in  a straiahtforward 
manner  by  partial  inteoration  with  the 
restriction  that  the  source  strenoth 
vanishes  at  the  body  ends. 

The  outer  and  inner  solutions  are  matched 
by  equating  (8) , (12)  and  (22).  The  ; 

terms  orooortional  to  log  r match 

automatically,  since  the  same  source  < 

strenqth  has  been  anticioatod  ^'or  the  ; 

two  solutions.  The  rem.aininq  terms,  of  ' 

order  one,  determine  the  "constant"  f.  (x) 
in  the  inner  solution  as  the  sum  of  the 
last  two  ‘■arms  in  (22a)  . It  follows  that 
the  inner  solution  (8)  can  be  expressed 


fr  L,  (Kx-KC)  dC  , (22a) 

J r» 


in  the  alternative  forms 


or 


S a(x)  [ 


log  (Kr)  + Y + mi 


(C)  I,^(Xx-KC)  d?  (22b) 


Here  the  kernels  L are  defined  as 
n 


L^(u)  •=  -log  {21  ul)  sgn(u),  (23a) 


(u)  ^ rHo(!ul)+  Yg(lul)+2i  J^duh]  , 


(23b) 


r f I 'il  1 

Lj  (u)  =|-log  2lu[-Y  -“i  Li(u')du'| 


J, 


' o J 

sgn(u)  , (23c) 


cind,  in  the  usual  notation,  H , Y and 

o o 


denote  the  Struve  and  Bessel  functions 


of  order  zero.  The  last  term  in  (22a) 
vanishes  as  K-0  , leaving  the  result 
for  a slender  body  in  an  infinite  fluid 
with  the  kernel  (23a)  governing  longitud- 
inal interactions  along  the  body.  The 
last  term  in  (22a)  can  therefore  be  inter- 
preted as  the  free-surface  correction,  and 
the  corresponding  kernel  (23b)  can  be 
derived  in  a direct  manner  by  subtracting 
(15)  from  (14),  followino  the  approach 


Ak  ' . 


= (y,z;x)-t-|-  j a(C)L,  (Kx-KC)dC  , 

A •'  0 


(24a) 


^ (2D)  (y,j,.^)^l_  j (x)  ^loq(K) +Yirij 

1 f' 

Jo 


+ =irr  1 a‘ (OLj  (Kx-KC)d?  . (24b) 


(2D) 


The  potential  4 is  the  solution  of 

the  two-dimensional  Laplace  equation  (7) , 
the  rioid  free-surface  condition  (11) , 
and  the  bodv-boundary  condition  (2).  The 
arbitrary  constant  in  this  solution  is 
determined  uniquely  by  the  requirement 
that  (12)  holds  at  large  distances  from 
the  body  section,  with  an  error  o(l). 

The  source  strength  o is  determined  by 
continuity,  as  twice  the  product  of  the 
local  waterolane  beam  and  vertical  velocity. 
The  only  effects  of  the  free  surface  are  on 
the  second  term  of  (24a) . 


There  are  two  practical  obiections  to  the 
above  theorv.  One  is  that  waves  exist 
only  in  the  outer  region.  The  other  is 
that  the  leading-order  force  opposing  the 
motion  of  the  ship  is  hydrostatic, 
proportional  to  the  waterolane  area  and 
hence  0(e).  Bv  comparison,  the  hvdro- 
dvnamic  force  and  the  inertial  force  due 
to  the  mass  of  the  ship  are  0(e‘).  Thus, 
to  leading  order,  there  is  no  resonance 
See  errata.  Page  115. 
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of  the  pitch  and  heave  motions.  A simple 
dimensional  argument  can  be  used  following 
Newman  (1977),  to  show  that  the  resonant 
frequencies  are  of  order  (g/B) ''^'=0(£  . 

These  defects  in  the  long-wavelength 
slender-body  theory  can  be  overcome  by 
including  higher-order  terms  in  e , 
following  Newman  and  Tuck  (1964).  Alter- 
natively we  can  revise  our  assumption 
regarding  the  order  of  magnitude  of  the 
frequency,  and  treat  the  short-wavelenath 
problem  where  1/L=0(£)  or  K=0(£~’). 

Short  Wavelengths,  K=0(£~M 

With  K=0(e~M,  both  terms  in  the  free- 
surface  boundary  condition  (3)  are  of  the 
same  order  of  magnitude  in  the  inner 
region.  Thus  we  anticipate  significant 
wave  effects  in  the  inner  region  adjacent 
to  the  body.  The  same  is  true  in 
the  outer  region,  with  the  assumptions 
that  radiated  waves  are  present,  of 
wavenumber  K , and  that  the  gradient 
is  proportional  to  K . Thus,  in  the 
short-wavelength  regime,  the  free- 
surface  condition  (3)  applies  in 
both  regions. 

Since  the  boundary-value  problem  for  the 
outer  solution  is  unchanged  from  the  long 
wavelength  case,  the  solution  can  be  con- 
structed in  an  identical  manner.  Re- 
stricting our  attention  to  vertical 
motions  of  ships  with  symmetry  about  th'=' 
plane  y=0,  the  outer  solution  is  given  by 
a distribution  (16)  of  three-dimensional 
sources  (14),  and  the  Fourier  transform 
of  the  outer  solution  is  given  by  (18). 
Only  the  transform  G*  must  be  reevaluated, 
with  K=0(£“M.  (Throughout  we  assume  that 
k=0(l),  or  that  longitudinal  variations 
of  the  quantities  transformed  are  qoverned 
by  the  ship's  length  as  opposed  to  the 
wavenumber.  This  assurption  is  not  valid 
in  the  far  field,  v/here  radiated  v/aves  of 
large  wavenuml^or  K prooacrate  in  all 
radial  directions.  However  the  Fourier 
transformed  quantities  are  to  be  used  onlv 
in  the  m.atchinu  reqion  where  this  assump- 
tion is  more  anorooria-r.  \ 


For  K>>1  , the  integral  (20)  for  C* 
can  be  approximated  by  deformino  the  con- 
tour of  integration,  following  Ogilvie 
and  Tuck  (1969)  and  Faltinsen  (1971). 

With  branch  cuts  of  the  scuare-root 
function  on  the  imaginary  axis  between 
±i<  and  infinity,  the  contour  of  inte- 
gration can  be  deformed  to  the  upper 
(lower)  half  plane  for  positive  (necativ'e) 
values  of  V . The  resulting  integral  is 
0(1/K),  and  the  dominant  contribution  is 
from  the  residue.  (Recall  from  (6)  that 
the  imaoinarv  part  of  K is  a small 
negative  quantity,  and  thus  the  poles 
of  (20)  are  adjacent  to  the  real  axis  in 
the  second  and  fourth  quadrants.)  TTie 
contribution  from  the  residue  gives 

G*  S -2mi  SXP  [K2-i!v'  (.K^-k^) 
d-k’/KM^'^- 

= -2mi  exp  [Kz-iKjyJ]  . (25) 

Using  the  last  form  o^  (25)  , which  is  V’alid 
for  K>^!k!  , the  inverse  transform  of  (18) 
is  given  in  the  Form 

♦ = — ic(x)exp  (Kz-iKlyl)  . (26) 

mhis  is  the  appropriate  inner  expansion 
the  outer  solution  ^or  the  short-wavelenoth 
case . 

The  inner  solution  is  governed  by  the  two- 
dimensional  Laplace  eouation  (7),  the  hody- 
boundarv  condition  (2) , and  the  free- 
surface  condition  (3).  Since  the  matching 
requirement  (26)  is  mathematically  identical 
to  a tv/o-dimensional  radiation  condition, 
the  desired  inner  solution  corresponds  to 
the  forced  motion  of  the  profile  C , with 
the  flow  constrained  to  the  plane 
x=constant.  This  problem  has  been  studied 
extensively,  for  a variety  of  body  pro- 
files. A compendium  o"  results  is  oiven 
bv  uucts  (lOf'8).  '-'ariou.-  computer  programs 
now  e.xist  for  solving  t.his  problem  wit.'i 
arhitrarv  body  profile  C,  as  ct.  scribed  bv 
Chapman  (1  ''77)  . 

Kith  the  inner  solution  specified  in  this 


i 
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manner,  the  problem  is  completely  solved. 
Since  the  frce-surfaco  condition  (3)  docs 
not  allow  an  arbitrary  constant  to  be 
added  to  the  inner  solution,  no  interac- 
tion exists  between  different  sections  of 
the  ship.  Thus  we  recover  the  strip 
theory  of  ship  motions,  with  the  flow  at 
each  section  of  the  ship  independent  of 
the  shape  or  motion  at  adjacent  sections. 

As  noted  in  the  Introduction,  the  last 
conclusion  is  cause  for  mixed  feelings. 

On  the  positive  side,  it  provides  a 
rational  basis  for  the  strip  theory, 
which  originally  was  introduced  to  ship 
hydrodynamics  in  a heuristic  manner.  On 
the  other  hand,  we  began  with  a fully 
three-dimensional  problem  and  it  is 
disappointing  to  find  that  in  the  final 
result  three-dimensional  effects  are 
absent  from  the  inner  solution.  More- 
over, if  forward  speed  is  included,  as 
in  Ogilvie  and  Tuck  (1969) , the  leading- 
order  solution  in  the  inner  region,  does • 
not  depend  on  the  ship's  forward  speed. 

As  in  the  analogous  situation  for  the 
long-wavelength  theory,  one  may  over- 
come the  deficiencies  of  the  simple  strip 
theory  by  including  higher-order  terms 
in  e . In  effect  this  is  done  in  all 
strip  theories  of  ship  motion  where  for- 
ward-speed effects  are  included.  Ogilvie 
and  Tuck  (1969)  derive  a consistent 
higher-order  t.heory,  including  terms  of 
relative  order  but  neglecting  0(e). 

Included  in  the  neglected  terms  are  con- 
tributions to  the  force  and  mom.ent  pro- 
portional to  the  square  of  the  forward 
speed,  which  are  assumed  generally  to  be 
of  practical  importance. 

This  situation  is  not  entirely  satisfactory. 
As  has  been  stated  aptly  by  Ogilvie  (1976) 
in  the  context  of  wave-resistance  theory, 

"we  are  more  likely  to  make  practical 
progress  by  developing  a better  linear 
theory  than  by  devising  second-order 
corrections  to  an  inadequate  linear 
theory" . 


A ComDOjite  (Approximation 
An  "interoolation"  between  the  two 
theories  outlined  above  has  been  con- 
structed by  Maruo  (1970) , in  a similar 
manner  to  the  additive  comoosition  of 
matched  asymototic  expansions.*  First 
we  note  that  the  inner  solutions  for 
K=0(1)  and  K=0(e~')  possess  a common 
"overlap"  expansion.  Thus  the  formal 
approximation  of  (24)  for  K is 

equal  to  the  approximation  of  the  strip- 
theory  potential  for  snail  K . The 
former  result  is  obtained  simply  by  noting" 
that  the  kernel  Lj  , defined  by  (23c) , 
vanishes  for  larae  values  of  its  argument. 
Thus  the  larae-K  approximation  of  (24) 
is  given  by 

'‘’(K=0(l))='*^^°l^a  (X)  [loa(K)+Y+^i]. 

(27) 

The  small-K  approximation  of  the  strip- 
theory  potential  can  be  obtained  by  ex- 
pressing the  solution  as  a source  plus 
wave-free  potentials,  following  Ursell 
(1968),  and  noting  that  the  second  term 
of  (27)  is  the  low-frequency  limit  of  the 
source  potential,  minus  the  loa(r)  term 
in  p (20)  ^ 


Following  the  rules  for  additive  composi- 
tion, we  add  the  two  inner  solutions  and 
subtract  their  common  limit  (27).  It 
follows 

4 = 4> 

s 

<7(x)  [log  (K)  + Y + TTi] 

“ *strip^  |^(r)L:(Kx-KC)dC  . (28) 

This  is  a composite  inner  solution,  with 
the  property  that  it  contains  the  K=«0(1) 
and  K=0(c”')  results  as  limitina  cases. 
However  (28)  does  not  satisfy  the  inner 
boundary-value  problem  for  all  values  of 
K , since  the  interaction  terms  added  to 

See  errata,  Page  115. 

* Computations  by  Maruo  are  reported  in 
discussion  of  Chang  (1977).  A similar 
approach  is  described  by. Tuck  (1966). 


that 
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I a(S)Lj  (Kx-KS)dC 
•'  o 
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the  strip-theory  potential 

violate  the  free-surface  condition  (3) . 

In  the  interpolation  theory  of  Maruo 
(1970),  the  interaction  terms  of  (28) 
are  multiplied  by  a higher-order  factor 
(l+Kz).  As  a result  the  free-surface 
condition  (3)  is  satisfied.  However,  the 
uniformity  of  this  correction  is  doubtful 
for  large  K,  and  a troublesome  term  which 
results  is  ignored  in  Maruo *s  computation 
of  the  hydrodyncunic  pressure  force. 

A Unified  Solution 

To  overcome  the  deficiencies  of  the  long- 
and  short-wavelength  theories,  it  is 
logical  to  seek  a "unified"  theory  which 
embraces  both  as  special  cases  but  is  valid 
throughout  the  wavenumber  regime 
0(1)<  K ? 0(e”*).  To  ensure  that  all 
leading-order  ingredients  of  the  two 
limiting  problems  are  included,  the 
guiding  priciple  is  that  boundary-value 
problems  are  formulated  in  the  inner  and 
outer  regions  such  that  for  all  relevant 
values  of  the  uavenurber  no  terms  of 
leading  order  in  the  slenderness  parameter 
are  neglected. 

Proceeding  on  this  basis  for  the  ship- 
motion  problem,  the  outer  problem  is  un- 
changed, but  in  the  inner  region  both 
terms  of  the  free-surface  condition  (3) 
are  retained  in  the  "inconsistent"  manner 
suggested  by  Ogilvie  (1974).  The  strip- 
theory  potential  of  the  short-wavelength 
regime  is  a particular  sol'tion  of  the 
resulting  inner  problem.  However  for 
arbitrary  wavelengths  it  is  not  possible 
to  match  this  solution  with  the  outer 
potential.  Indeed,  for  long  wavelengths, 
an  arbitrary  additive  constant  must  be 
allov/ed  in  the  inner  solution,  as  in  (8)  . 
Since  the  free-surface  condition  (3)  does 
not  permit  such  a constant  to  be  added  to 
the  inner  solution,  a nontrivial  homo- 
geneous solution  is  required. 

The  only  solutions  of  the  homovTcnoous 
problem  in  the  inner  rc.cicn  are  the 
velocity  potential.;  of  the  diffraction 
problem,  for  scattvrinq  of  incident  waves 
m.oving  past  the  fi-xcU  body  in  the  oositive. 


and  negative  directions.  The  appropriate 
solution  symmetrical  about  y*0  is  com- 
posed of  two  equal  and  opposite  incident 
waves,  or  an  "incident  standing  wave". 

This  homogeneous  solution  violates  the 
two-dimensional  radiation  condition  but, 
from  (26) , that  objection  applies  only 
for  short  wavelengths.  As  a result, 
we  anticioate  that  the  contribution  to 
the  inner  ootential  from  the  homoaeneous 
solution  must  vanish  in  the  short-wave- 
lenoth  reaime.  More  generally,  a non- 
zero contribution  from  the  homogeneous 
solution  is  both  necessary  and  aDorooriate 
to  match  with  the  outer  solution. 

Thus,  the  unified  slender-body  theory  for 
ship  motions  in  calm  water  is  based  on 
retaining  the  complete  free-surface 
condition  (3)  in  both  the  inner  and  outer 
regions.  The  inner  problem  is  simolified 
principally  by  imposing  the  two-dimensional 
Laplace  equation  (7) , which  is  valid  for 
all  wavenumbers  provided  waves  in  the 
inner  region  are  parallel  to  the  x-axis. 

The  inner  solution  is  assumed  to  consist 
of  a linear  combination  of  the  two- 
dimensional  strip-theory  ootential, 
correspondina  to  forced  motion  of  the 
body,  and  of  the  two-dimensional  hom.oaene- 
ous  solution  for  diffraction  of  incident 
standing  waves  by  the  fixed  bodv.  For 
the  outer  solution  a longitudinal  distri- 
bution (16)  of  three-dimensional  sources 
(14)  is  assumed.  The  resultina  inner 
and  outer  solutions  can  be  matched,  in 
a suitable  overlap  region  B<<;y<<L 

This  unified  theory  of  ship  motions  can 
be  compared  to  the  modified  strio  theory 
of  Grim  (1960).  Grim  solves  the  two- 
dimensional  problem  with  acorooriate 
sinaularities , and  assumes  a corresnondi.-.c 
distribution  of  three-dimensional  sinau- 
larities with  the  same  coefficients.  To 
correct  the  resultina  error  in  the  body 
bcundarv  condition,  duo  to  the  change  from 
uv/o  to  th'-oo  dir’en.sions  , a ' t r 1 

diffraction  tolution  is ’ introduced . I!ero 
the  incitient  waves  are  associated  with 
radiation  along  the  hull  from  t!ie  three- 
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dimensional  wave  source.  Grim  anticipates 
that  this  "longitudinal  diffraction 
problem"  can  not  be  realized  physically, 
as  has  been  confirmed  mathematically  by 
Ursell  (1975). 


The  inner  region  is  defined  by  the  scale 
of  the  maximum  body  radius,  r=0(c).  Here 
a coordinate  stretchinq  argument  can  be 
used  to  reduce  (29)  to  its  two-dimensional 
form 


Aside  from  the  lack  of  a systematic  form- 
alism based  on  asymptotic  approximations, 
and  the  introduction  of  additional 
assumptions  by  Grim  to  simplify  his  comp- 
utations, the  principal  distinction  with 
the  present  unified  theory  is  that  we 

supplement  our  inner  solution  with  a 
/ 

diffraction  potential  for  transverse. 
incident  waves.  The  resulting  solution 
is  strictly  two-dimensional,  and  well 
defined  both  mathematically  and  physically. 
The  expected  longitudinal  interactions 
along  the  hull  arise  indirectly  from 
matching  with  an  outer  solution. 

The  analysis  of  ship  motions  based  on  the 
unified  slender-body  theory  is  now  in 
progress,  both  with  and  without  forward 
speed.  Rather  than  pursuing  either  study 
here  in  detail,  I shall  describe  instead 
a simpler  but  analogous  problem  where  the 
same  approach  can  be  illustrated  more 
exp] icitly . 

3.  AN  analogous  ACOUSTIC  PftOBLE.M 
Let  us  consider  a slender  axisymmetric 
body,  of  radius  r»r^(x),  which  radiates 
acoustic  waves  in  the  surrounding  medium. 
Nith  the  usual  assumptions  of  linearized 
acoustics,  the  Laplace  equation  (1)  is 
replaced  by  the  reduced  wave  equation 


= 0 , (31) 

yy  zz 


with  the  error  a factor  1+0 (e*)  as  in 
(7)  . The  most  general  axisymmetric 
solution  of  (31)  is 

()  = a (X)  H^^^’(Kr)  + a,  (x)  H^<^’(Kr)  , 

1 o z o 

(32) 


where  a,  and  a,  are  coefficients 

* (12) 
which  may  depend  on  x , and  ' 

denote  the  Hankel  functions  J^iiY^ 

o o 

With  the  assumed  time  dependence 

e^“^,  H and  H represent  in- 

o o 

coming  and  outgoinq  waves,  respectively, 
for  large  Kr. 

The  outer  solution  is  a longitudinal 
distribution  of  acoustic  sources  along 
the  body  axis.  Thus,  as  in  (16) 


G(x-5,r)  dS 


where 


G(x,r)= 


1/5 

^-iK(x*+rM 
e 

(x*+ 


(33) 


(34) 


(V*  + K^)  - 0 , (29) 

where  K=»ojL/c  is  the  nondiraensional  wave- 
number  and  c is  the  speed  of  sound.  The 
body  boundary  condition  (2)  is  unchanged, 
and  the  radiation  condition  of  outgoing 
waves  at  infinity  can  be  replaced  by  the 
complex  wavenumber  (6). 

For  simplicity  we  assume  axisymmetric 
forced  motions  of  the  body,  with  normal 
velocity  V(x).  The  leading-order  boundary 
condition  on  the  body  takes  the  form 

■ -V(x),  on  r » r^(x).  (30) 


Fourier  transforms  can  be  used  as  in 
Section  2,  to  derive  the  inner  expansion 
of  the  outer  solution,  for  purposes  of 
matchinq  with  (32) . Here  we  follow  a 
modified  approach,  matching  the  trans- 
forms of  both  the  inner  and  outer  solutions. 
For  this  purpose  the  transform  of  the 
inner  solution  (32)  is 

$*=«  a,*  H (Kr)  + a,*  H (Kr)  . (35) 

o ^ o 

Using  the  convolution  theorem  (18) , 
evaluating  G*  from  (34),  and  relating 
this  to  the  intooral  representation 


of  the  Hankel  function,  the  transform  of 
the  outer  solution  is 


f2  K^(r/k=-K*) 


-iri  H (r/K=^0 
o 


f for  {IIc|<K}. 


(36) 


Matching  now  can  be  performed  in  the 
transformed  domain,  by  equating  (35)  and 
(36)  for  a suitable  ranae  of  r in  the 


overlap  domain  r <<r<<  1.  As  in  Section 
o 


2,  we  assume  that  k»0(l),  and  e<<kr<<l 
in  the  matching  region. 

Short  Wavelengths,  K»0(e~*) 

For  short  waves,  where  K^>lk[  , the 
second  form  of  (36)  applies.  After  ex- 
panding for  small  values  of  [kj/K  , 
the  transformed  outer  solution  is 


**  - -h  [»c 


(2) 


(Kr) 


+ j (k*rVKr)Hi  (Kr)j  , 


(37) 


where  the  error  is  a factor  l+0(k'’r*/K* ) . 
Equations  (35)  and  (37)  can  be  matched  to 
leading  order  with  the  results 


a,*=0 


* * i ia* 


(38) 

(39) 


and  the  error  is  a factor  l+0(k'r/K). 
Transforming  (38-39)  determines  the 
coefficients 


a,  - 4 la 


(40) 

(41) 


Thus,  the  inner  solution  for  short  wave- 
lengths is  given  to  leadino  order  by 


;=a.  (X)  fi  (Kr)=  - 


V(x) 


111 


(2)  - (Kr)  . 

KH,  (2)(Kr  ) ° 


where  the  last  equality  follows  from 
the  boundary  condition  (30).  This  is 
the  "strip  theory"  result,  with  the 
inner  solution  identical  to  the  two- 
dimensional  radiation  from  a circular 
cylinder.  Since  ai»0  , the  two- 
dimensional  radiation  condition  is 
satisfied  in  an  analoqous  manner  to 
the  strio  theory  of  ship  motions. 

Long  Wavelengths,  K°0(1) 

For  (k,K)«0(l),  and  small  r,  the  Hankel 
functions  in  (35)  and  (36)  can  be  approxi- 
mated from  the  leading  terms  in  their 
infinite-series  reoresentations.  .^'atching 
the  two  transformed  solutions  after  this 
is  done  gives  the  result 


(ai*+aj*) (1-^  K*r*)  + | i(a,*-a*)  tl-jK-r=] 
|^log(i  Kr)  + Y + j K*r*| 


— o*  ri+-(k^-KMr* 
1-4 


271 


VJ  Kr)  + Y + j log(IkVK'-lI 


[loo  4 

^ (k--K>)r-  .{Si/,}],  !orf|k!;K}  . 


(43) 


Here  y * 0.577 is  Euler's  constant 

and  the  error  in  (43)  is  a factor 

1+0 ( (kr) " , (Kr ) * ) . To  leading  order  for 

small  r , it  follows  that 


1 f log  (k’/K’-l)  ] 

a, * + aj*  = ^ [ 

* llogd-kVkM  + TTi]  , 


a * = - ± ic* 
2 4 


(44) 

(45) 


where  the  error  is  a factor 
1+0  ( C-cr)  • , (Kr)  ")  . 


The  inverse  uransforn  of  (45)  is  simply 


(42) 


^T^BSPMPIvssr^MHPW — liMi  I 
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a,  - aj  = -j  io  . (46) 

The  inverse  transform  of  the  factor  in 
braces  in  (44)  can  be  related  to  the 
exponential  integral,  after  a contour 
integration  in  the  complex  )c-Dlane  with 
branch  cuts  on  the  real  axis  for  |lcI>K 
Another  application  of  the  convolution 
theorem  then  gives  the  transform  of  (44), 

a,  + a,  = J io  + ^ |‘c' (C)L(Kx-KC) )dC  . 

(47) 

Here  the  )cernel  L is  defined  by 

L(u)  = Ej(ilul)  sgn  (u)  , (48) 

and 

E.  (z)  =«  [ dt  (49) 

•'  z t 

is  the  exponential  integral 

After  solving  (46)  and  (47)  for  the 
coefficients  a and  a,,  and  substitut- 
ing  in  (32) , the  inner  solution  can  be 
expressed  in  the  form 


(52)  in  (33) . This  asocct  of  the  lona- 
wavelenath  oroblem  has  been  studied  bv 
Pond  (1966)  , who  deduces  the  source 
strenoth  (52)  in  a alobal  manner  without 
consideration  of  the  inner  solution  for 
the  velocity  potential.  Chertock  (1975) 
has  presented  a more  aeneral  s lender-body 
theory  'or  the  lona-wave length  case, 
starting  with  an  integral  equation  for 
the  exact  three-dimensional  solution  and 
using  physical  arguments  to  approximate 
the  three-dimensional  singularity 
distribution.  No  attempt  is  made  to 
reduce  the  inner  solution  to  the  form 
(51). 

For  ver'y  long  wavelengths,  K<<1  , the 
kernel  (48)  can  be  approximated  from  the 
infinite-series  expansion  of  (49) . In 
the  limit  K-^0  , (51)  tends  to  the 
classical  inner  solution  for  an  unbounded 
incompressible  fluid, 

+ |7[[cT'(C)  log(2lx-El)sgn(f-x)df  . 

(53) 


|i 

I 


p ia  (X)  (Kr) 


This  is  consistent  with  the  corresponding 
limit  of  (22)  . 


1 n 

+ ^ Jo(Kr)j^a' (5)L(Fj<-KS))d5.  (50) 

Finally,  if  the  Bessel  functions  in  (50) 
are  expanded  for  sm.all  Kr,  it  follows  that 

d “ ^ a(x)j^log(2-  Kr)  + Y + jij 

+ J7|V(5)  I.(Kx-KC))d?  . (51) 

At  this  stage  the  boundary  condition  (30) 
can  be  used  to  determine  the  source 
strength 


A Composite  Approximation 
The  inner  solutions  (42)  and  (51)  are 
valid  in  the  mutually  exclusive  reaimes 
K“0(e“  ) and  K*0(1) , respectively. 

These  solutions  possess  a common  approxi- 
mation in  the  "overlap  reoion"  l<<K<<e~’  , 

- r^V  [log(j  Kr)  + Y + 7 i]  . (54) 

which  follows  by  assuming  Kr<<l  or 
K>>1  in  (42)  or  (51),  respectively.  Thus 
a composite  solution  can  be  derived  by 
addition  of  (42)  and  (51),  and  subtraction 
of  (54),  in  the  form 


0 (x)*  -2!Tr^  (x)V(x)  , 


(52) 


and  the  inner  solution  is  completed. 

The  outer  solution  and  far-field  radiation 
patterns  can  be  derived  by  substituting 


V (x) 


(2) 


K H.  ‘‘^Kr  ) 
o 


(Kr) 


h f 

) 0 


) L(Kx-K")d£: 


(55) 
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with  a defined  by  (52).  The  analooous 
formula  for  the  shio-motion  oroblem  is 
(28).  Once  again  the  composite  solution 
is  correct  in  the  two  limiting  regimes 
where  K»0(1)  or  0(£”’)  , but  the 
validity  is  not  uniform  throughout  the 
intermediate  ranoe  of  wavenumbers  since 
the  last  term  in  (55)  violates  the 
wave  equation. 


A Unified  Solution 

The  deficiency  in  (55)  can  be  overcome 
simply  by  observing. that  the  results 
of  our  derivation  for  the  long-wavelength 
case,  up  to  and  including  (50), are  con- 
sistent with  the  corresponding  short- 
wavelength  results.  Thus,  for  |lcI/K<<l  , 
the  Fourier  transforms  (44-45)  yield 
(38-39).  Similarly,  since  the  integral 
in  (47)  tends  to  zero  for  K>>1,  (46-47) 
yield  the  coefficients  (40-41)  and  the 
inner  solution  (50)  is  dominated  by  the 
first  term,  proportional  to  H^^2)  ^ 
agreement  with  (42). 

From  this  standpoint  the  inner  solution 
(50)  apparently  is  valid  for  all  wave- 
numbers  KfO(£“*)  , provided  the  source 
strength  is  determined  in  accordance 
with  the  boundary  condition  (30).  This 
gives  an  integro-dif ferential  equation 
for  the  source  strength, 

ic (X)  H,  (Kr^(x) ) 

^ I (^)L(Kx-Kf)  )d-;  = i V(x)  . 

(56) 

To  confirm  that  (50)  is  the  desired 
"unified"  solution,  we  note  that  this 
velocity  potential  satisfies  the  two- 
dimensio.-.al  wave  equation  (31)  and  the 
boundary  condition  (30}  in  the  inner 
field;  the  corresponding  outer  solution 
(33)  satisfies  the  three-dimensional 
wave  equ.Jticn  (29)  and  the  aonront-i  at  o 


radiation  condition.  The  only  question 
which  remains  is  to  confirm  that  (33)  and 
(50)  can  be  matched  for  all  oossible 
wavenumbers,  in  some  appropriate  overlap 
region  close  to  the  body  in  the  outer 
problem  and  far  from  the  body  in  the 
inner  problem. 

The  overlap  reaion  is  expected  to  exist 
for  a range  of  r such  that  e<<r<<l, 
or  when  r»0(c°‘),  0<a<l  . Since  the 

Four ier-trans form  parameter  lc“0(l), 
ltr*o(l).  However  the  value  of  Kr  in 
the  overlap  domain  is  not  restricted  to 
this  range.  Xn  particular,  when  the 
body  is  radiating  short  waves,  the  match- 
ing must  be  carried  out  in  the  .far-field 
of  the  inner  solution,  where  Kr>>l 
Since  the  matching  condition  (43)  neglects 
a factor  l+0(Kr)'‘,  and  the  unified  inner 
solution  (50)  derives  from  this  condition, 
the  possibility  exists  that  (43)  is  not 
valid  throughout  the  range  of  all  wave- 
numbers  K<0(e“’),  in  spite  of  the  fact 
that  it  yields  correct  results  in  the 
limit  K»0(c“’).  This  ouestion  can  be 
resolved  by  noting  that  the  error  in  (43) 
is  a factor  l+0(Kr)*,  whereas  the  error 
in  the  short-wavelength  results  (38-39) 
is  a factor  l+0(lc^r/K)  . Requiring  these 
errors  to  be  of  equal  maanitude  gives  a 
condition  for  transition  between  the 
short-  and  long-wavelenath  matchino 
relations , 

k^r/K  = (Kr) ‘ . (57) 

If  r=0(z'^)  , the  transitional  wavenumber 
is 

V = --sa/.  (58) 

K = e.  , 

and  the  error  in  matchina  is  a factor 

l+0(c'‘^'^’)  . Thus  the  matchino  error  will 

be  0(£)  if  the  overlap  rcoion  is  defined 

by  r=0(c’^’)  at  the  transitional  wavenumber 
-V, 

K = c .A  more  comolete  analysis  is 
necessnrv  to  prove  that  the  full  solution 
is  accurate  to  0(r)  with  this  oarticular 
choice  o*  the  overlao  reoicn.  Neverthe- 
less it  is  apparent  that  the  relative 
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error  in  the  unified  inner  solution  (30) 
is  o(l)  for  all  wavenumbers. 

The  unified  solution  (30)  can  be  compared 
with  the  composite  approximation  (35). 

The  latter  result  is  equivalent  to  re- 
placing the  Bessel  function  J^(Kr)  by 
its  limit  of  unity,  in  the  second  term 
of  (50) , and  determining  the  source 
strength  o in  each  term  of  (50)  from 
(42)  and  (52)  respectively. 
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Errata 

The  term  in  (22a)  involving  the  kernel  has  been  deleted  from 

(24a)  and  (28) . To  correct  this  error  the  second  line  of  (22a) 
should  be  inserted  before  the  plus  sign  in  (24a)  , and  before  the 
first  plus  sign  in  (28) . 


